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Using iso-spin symmetry the quark-connected and disconnected contributions to the hadronic 
vacuum polarisation in a theory with Nf = 2 flavours can be described as independent correlation 
functions, respectively. We show how this allows to use twisted boundary conditions for the 
connected contribution in order to improve the q 2 -resolution in lattice QCD. Furthermore we 
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chiral perturbation theory. We discuss extensions to theories with more than 2 dynamical flavours. 
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Table 1: Some recent results for the leading hadronic contribution to the muon anomalous moment, a™ 

1. Introduction 

As recently summarized in ^ the current experimental value for the full muon anomalous 
magnetic moment is = 11659208(6.3) x 10~ 10 which has to be compared to the Standard 
Model prediction of = 11659179(6.5) x 10~ 10 . Although there are discussions about how to 
correctly estimate the size of the systematic errors there persists a tension between experiment 
and theory which amounts to around three standard deviations (3.2a in the comparison quoted 
above). The muon anomalous moment is a remarkable observable in that all three sectors of the 
SM contribute considerably. While perturbation theory is employed in order to predict QED and 
weak contributions, the leading and next-to-leading hadronic contributions are non-perturbative 
effects. Lacking reliable and precise theory computations from first principles the current SM pre- 
dictions for the leading hadronic contribution ajf d are derived from experimental measurements 
of e + e~ -annihilation into hadrons leading to a world average ^ for the hadronic contribution of 
fl had = 690.3(5.3) x 10- 10 . 

Table [j] summarizes previous attempts to compute a^ A in lattice QCD in comparison with 
the determination from e + e~ annihilation. Without going into the technical details of the various 
attempts to compute a h ^ d on the lattice it is obvious that currently the purely theoretical predictions 
can't match the level of precision which can be reached via the experimental determination. But 
yet, an independent confirmation of the result from e + e~ -annihilation matching it in precision and 
by means of a SM calculation is clearly desireable. 

The leading hadronic contribution to the muon anomalous moment is the convolution integral 

DO 

a jad = 4n 2 ^y j dK 2 f{K2) ( n(^) _ n(0 )) , (i.i) 







where K is the Euclidean momentum and the function f{K ) in (1. 1) diverges for K 2 — > (see e.g. 



the discussion in |Q]). H(K 2 ) is the vacuum polarization which for a theory with Nf quark flavours 
is defined through 

n% f \q) = ( q , qv - g , v qW) = if d^^j^ EM \x)ji N f- EM \0)} , (1.2) 

where ja f,E is the corresponding electromagnetic current. 

Apart from the usual sources of systematic errors in simulations of lattice QCD (unphysically 
heavy quark mass, finite volume, cut-off effects) we identify three sources of systematic errors that 
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have to be brought under control in order to be able to make precise predictions for a^ A from first 
principles: 

1) since f(K 2 ) in ( |L~I| ) diverges for K — > it is crucial to have a good momentum resolution 
for IT for small values of K. Given that the smallest Fourier mode in eqn. ( |1.2| ) on the lattice 



for typical simulation parameters is jj- 250MeV (L being the spatial extent of the lattice 



2L 

2\ 



and usually the time-extent is chosen as T = 2L) one currently extrapolates results for Tl(q 
from the region where lattice data is available using some phenomenological ansatz 



2) carrying out the Wick contraction in (1.2) reveals contributions from quark-disconnected 



diagrams which are notoriously hard to compute in lattice QCD 

3) Aubin and Blum ^ found strong signs for vector dominance in their lattice data which has 
to be modelled in order to allow for a reliable description of the data within an effective 
theory frame work. 

Here we suggest a computational strategy that should allow to systematically reduce the uncertain- 
ties due to 1) and 2). Concerning 3) we will follow [Q] as long as the simulated quark masses are 
unphysically large. They used a model where a vector particle is added to the chiral Lagrangian as 
an additional degree of freedom in order to better describe their lattice data for Yl{q 2 ). This work is 
part of a major effort by the Mainz group aiming at a precision computation of the leading hadronic 
contribution to the muon anomalous moment (see also Hartmut Wittig's talk [|6|]). 

2. New strategy for Nf = 2 

Partially twisted boundary conditions [J7|, ||[ ^ 11, |l2|, 13, 14, 15] have now been used 
successfully to improve computations of observables of processes that depend on the hadron mo- 
mentum. However, as already stated in [Q] the net effect of the twist in flavour neutral hadrons, 
like e.g. the n is zero since the twist of the quark and the anti-quark cancel. The situation is sim- 
ilar here - the electromagnetic current is flavour diagonal, hence a naive application of (partially) 
twisted boundary conditions is ruled out. 



The correlator entering eqn. ( |1.2| ) for Nf = 2 is 



where we make the flavor content of the quark bilinear currents explicit and where we have fac- 
torized out the electromagnetic charges of the quarks in units of e on the r.h.s.. The first and the 
last correlator on the r.h.s. receive contributions from both a quark-connected and a disconnected 
diagram and the second and third correlator are disconnected, only. After carrying out the Wick 
contractions and using iso-spin symmetry (m q = m u = mi), the two point function one has to com- 
pute in a lattice simulation is 



C { -f M \q) = ^ x {(^r{q{x)Y v q{x)q{0)y il q{0)} + 

x (2 2) 

±Tr{y w (x)?(x) }Tr{ y^(0)?(0) })} , 
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which again using iso-spin can be written in the equivalent form 

cf v EM \q) = (Tr{u(x)r v d(x)d(0)^u(0)} ') + 



„(2,con) , , 

. _ r ... ,-^\^ (2-3) 

9 



Tr { j v d(x)d{x) } Tr {y M M(0)w(0) } \ } . 



=4 2 / isc) ( 9 ) 

By using iso-spin we are able to express the vacuum polarization in terms of two correlation func- 
tions C^v° n ^ (q) and C^ v dlsc ^ (q) with their individual continuum and infinite volume limits. The 
crucial point to note is that in this way the coupling of the photon to the quark in the connected 
piece has been replaced by a coupling to a flavour non-diagonal current and partially twisted bound- 
ary conditions can be applied. We note that a similar trick has been used for computations of the 



electromagnetic pion form factor in [llq, 17] 



We are still left with the question of how to treat the disconnected part. To this end we resort 
to the description of the vacuum polarization in SU (2) chiral perturbation theory at NLO. In SU (2) 
the electromagnetic current receives the following singlet and non-singlet contributions, 

(2 4) 

if = 2^(01 +iC2)7iiV> 4" = -to) 7^ V, 

where \(f T = (u,d), (T; for i = 1 , 2, 3 are the Pauli matrices and Oq = y / l/2l2 X 2- The effective theory 
that describes the Nf = 2 vacuum polarization is SU (2) chiral perturbation theory as formulated in 
[ |T8| ] but including at 0(p 4 ) terms with non-vanishing flavour-trace, 

JSfW = 2l 5 (U% v U& v ) +Ah 2 {R iiV R^ +L^ V ) +4h s (R^ v +L^)(R^+Ln , (2.5) 

V V 

flavour off— diagonal flavour— diagonal 

where angular brackets indicate a flavour trace and the hat indicates that the latter has been sub- 
tracted (cf. Jl9| ] for further details on the notation). Note the additional low energy constant h s 
multiplying the flavour diagonal term. 

The diagrams contributing at NLO are illustrated in figure |] and our preliminary result is 

n$( q )=n$ 3 Hq) + ln$°\q) 

= -(^^v-^v9 2 )('4B2l(9^m2) + 2/5(A^o)+4/^ 2 ( i ^o) + 5/^s + 4^ ? ln(^|)) , 

where Hq is the renormalization scale and B2i(q 2 ,m 2 z ) is a loop integral (cf. [19]). The contribu- 
tions Tiny {a = 0, 1,2,3) are the effective theory descriptions of 2pt-correlators constructed of the 
currents c\j/y^/ v a a ^f, where c 6 C is a normalization. We add that n^ ^ (q 2 ) = n(° ! °) is momentum 
independent at NLO and only depends on the low energy constant h s . 

Alternatively it is possible to express the quark-connected and disconnected contributions in- 
dividually in the effective theory (the linear combination of which should of course again yield 
eqn. (2J5)). We computed the corresponding expressions as 

n (2,con) (9 2 )= 10^3,3)^2) 5 

(2.7) 

n (2,disc) (? 2 )= i(n(°.o)_ n ( 3 . 3 V)) . 
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W/^ ^^lA/W wwvwtwvvwv vwvw (J) www 
Figure 1: Diagrams in chiral perturbation theory that contribute to the vacuum polarization at NLO. 

Since n' ' is independent of q 2 , the quark-disconnected part turns out to have the same momen- 
tum dependence as the connected part up to a finite shift proportional to h s , 

n (2,disc) ^ = i n (o,o) _ i.n( 2 .™) (q 2 ) . (2 .8) 



Assuming the validity of chiral perturbation theory at NLO the knowledge of h s allows to fully 
predict the quark-disconnected contribution. 

The leading hadronic contribution to the muon anomalous moment is computed from 

tl(q 2 ) = l% 2 )-n(0). (2.9) 

One immediate observation is, that at NLO in chiral perturbation theory this expression is free 
of low energy constants since the loop integrals carry all the momentum dependence. A crucial 
consequence which to our knowledge hasn't been realized before is, that e.g. in the case of SU (2), 

falcon {ZAV) 

This is a remarkable result which tells us that at this order of the effective theory the disconnected 
part shifts the central value of the connected part by minus ten per cent. 



3. Conclusions and outlook 



To summarize, we have shown how partially twisted boundary conditions can be used to im- 
prove the momentum resolution in lattice computations of the hadronic vacuum polarization of a 
photon. It was shown how to analytically predict at NLO in chiral perturbation theory the contribu- 
tion of quark-disconnected diagrams to the leading hadronic contribution of the muon anomalous 
magnetic moment. It turns out that it's effect is to reduce the contribution of the connected part by 
10%. 

We are currently working on extending these arguments to the case of Nj ■ = 2 + 1 flavours. One 
complication there is that the iso-spin argument that allowed us to write the quark-connected and 
disconnected pieces as individual correlation functions in the case of Nf = 2 does not work straight 
forwardly since the strange quark doesn't have an iso-spin partner. Naively the ^^-contribution to 
the EM current can therefore not be written in terms of a flavor off-diagonal current. Our strategy 
is to extend the flavour group from SU (3) to SU (3 + 1| 1), i.e. to a graded flavour group with an 
additional partially quenched quark (which we call r-quark) that behaves like a mass-degenerate 
iso-partner of the s-quark. Then the quark-connected and disconnected pieces can be written in 
terms of individual correlation functions. 
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Since quark-disconnected diagrams are extremely difficult to compute on the lattice, we will 
investigate if our method to predict their contribution in the frame work of chiral perturbation 
theory can be applied to other phenomenologically interesting observables in QCD. 

Following [^] we are also computing the finite volume effects for the vacuum polarization 
including the effect of partial twisting. 

Another task that remains to be finished is the inclusion of vector degrees of freedom in order 
to be able to extrapolate the lattice data to physical quark masses and to assess the vector particle's 



effect on the stability of eqn. ( |2.10| ). 

As a joint effort the Mainz group has implemented correlation functions relevant for a lattice 
computation of a h ^ A in a C-code. First computations on gauge configurations of Nf = 2 non- 
perturbatively improved Wilson fermions which were generated as a collaborative effort within 
CLS (cf. the talks by Hartmut Wittig's ^ and Stefan Schafer [20] at this conference) have been 
carried out on the Mainz Wilson-Cluster. 

Acknowledgements: We warmly thank Dalibor Djukanovic, Jambul Gegelia, Chris Sachra- 
jda, Stefan Scherer and Hartmut Wittig for valuable discussions. 
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